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I.  INTRODUCTION 


Every  product  eventually  fails  and,  sooner  or  later,  is 
no  longer  usable  for  its  original  purpose.  The  rate  at  which 
this  deterioration  process  proceeds  is  called  the  failure  rate. 
One  is  interested  in  predicting  the  time  to  failure  for  a 
component  or  an  entire  system.  It  is  impossible  to  predict 
the  time  to  failure  precisely,  but  one  can  model  this  process, 
viewing  the  time  to  failure  as  a  random  variable  which  has 
a  probability  distribution  that  can  be  specified.  There 
exist  numerous  life  distributions  which  provide  realistic 
probabilistic  descriptions  of  the  behavior  of  components  or 
whole  systems. 

This  paper  is  concerned  only  with  the  exponential  distri¬ 
bution  which  plays  a  fundamental  role  in  reliability.  This 
distribution  has  many  properties  which  facilitate  the  modelling 
process,  and  so  encourage  its  application  in  many  fields. 

Where  not  directly  applicable,  it  is  often  used  as  a  bound. 

According  to  the  exponential  distribution,  a  component 
has  a  constant  failure  rate.  An  equivalent  property  is  that 
the  remaining  life  of  the  component  is  independent  of  its  age: 
the  "memoryless "  or  "no  aging"  or  "no  wear"  property.  This 
means  that  a  used  component  is  always  as  good  as  a  new  one. 

Convolving  life  distributions  for  components  to  get  the 
final  life  distribution  for  the  system  is  rather  difficult. 
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even  for  the  exponential  distribution.  For  an  increasing 
nunber  of  components  it's  no  longer  possible  to  compute  a 
closed  fora  life  distribution.  This  paper  will  show  how  to 
use  a  reliability  shorthand  notation  to  get  analytical  results 
for  small  basic  systems,  and  later  use  an  algorithm  to  com¬ 
pute  numerical  results  with  an  Apple  11+  microcomputer.  The 
Pascal  programs  are  provided  in  the  appendices. 


II.  RELIABILITY  SHORTHAND 


A.  NOTATION  AND  DEFINITIONS 

The  shorthand  notation  "EXP[A^]  +  EXP  [ A 2 ] "  denotes  the 
distribution  for  a  random  variable  T^+T2.  T^  has  an  exponen¬ 
tial  distribution  with  failure  rate  A^  and  T2  has  an  exponen¬ 
tial  distribution  with  failure  rate  X2«  T^  and  T2  are  assumed 
to  be  independent. 

In  the  following  T  always  denotes  a  random  life  length 

and 

F(t)  «  Probability  (T  >  t)  »  1  -  F(t) 

is  the  survival  function  as  a  function  of  time. 

The  convolution  of  life  distributions  corresponds  to 
the  summation  of  components  with  independent  random  lives 
to  get  the  life  distribution  of  a  system. 


III.  RELIABILITY  FOR  A  ONE  COMPONENT  SYSTEM  WITH  SPARES 


A.  GENERAL  DESCRIPTION 

There  are  two  ways  for  a  simple  system,  consisting  of 
one  active  component  and  one  spare,  to  survive  a  specified 
time  t.  Either  the  active  component  survives  time  t,  in 
which  case  the  spare  need  not  be  used,  or  the  active  com¬ 
ponent  fails  at  some  time  s  (0  ^  s  £  t)  and  the  spare  takes 
over  and  doesn't  fail  in  the  interval  from  s  to  t  and  thus 
completes  the  mission. 

Let  F^(t)  and  F2(t)  be  the  exponential  survival  functions 
for  component  one  and  two  respectively,  and  f^(t)  and  f2(t) 
their  respective  densities.  Then  the  survival  function  F(t) 
for  the  system  is  denoted  bys 

_  „  t 

F (t)  =  F. (t)  +  /  F, (t-s) f. (s)ds 

1  0  Z 

where  the  integrand  shows  the  probability  that  component  one 
fails  at  time  s  (f^(s))  and  component  two  takes  over  for 
at  least  the  length  of  time  t-s.  The  following  shows  some 
examples  how  to  get  the  survival  function  in  closed  form  for 
different  basic  convolutions. 
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B.  BASIC  SURVIVAL  FUNCTIONS 

1.  System  with  N  Identical  Components 

The  shorthand  notation  for  this  system  is: 

EXP [A]  +  EXP [A]  +  ...  +  EXP  [A] 


Whenever  an  active  component  fails  a  spare  takes  over,  as 
long  as  there  is  a  spare.  When  the  n-th  component  fails,  the 
whole  system  stops  functioning. 

For  n  =  2 : 

_  _  t 

F (t)  *  F.  (t)  +  /  F_  (t-s)  f .  (s)ds 

1  o  z 

*  e“Xt  +  f  e~X(t's)Ae"Xsds 
0 

-  (1  +  At)e"Xt,  t  >  0. 


The  density  function  for  this  survival  function  is: 

dF1(t) 


f  (t) 


St 


d[(l+At)e"XtJ 

- at - 


[ (1+At) (-Ae"Xt)  +  e_XtA] 


,  .  -At  , 2  -At  .  . -At, 
[-Ae  -  A  te  +  Ae  ] 


A2te”Xt 
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The  survival  function  for  the  two  component  system 
will  be  convolved  with  the  survival  function  of  the  next 
component  to  obtain  the  survival  function  for  the  three 
component  system. 

F(t)  =  (l+At)e“Xt  +  /fc  e"X  (t“s)  X2se"Asds 

0 

which  leads  to: 

2 

F(t)  *  (It  At  +  )e~Xt,  t  >  0. 

For  all  n: 

Continuing  in  this  manner  one  gets  the  general  solution 
for  the  survival  function  of  this  kind  of  system: 

n  4  1 

l  Ut)1  1 

F(t)  =  — - ryr—  ,  t  >  0. 

(i-l)e  At 

This  class  of  convolutions  always  gives  the  well-known 
Erlang [n, A]  survival  function. 

2.  System  with  N  Dissimilar  Components 

The  shorthand  notation  for  this  system  is  denoted 


EXP  [  A .  ]  +  EXP[A,J  + 


+  EXP [Al 


For  n  ■  2: 


The  survival  function  for  the  active  component  is 


Fx(t)  =  e 


-Xjt 


t  >  0. 


The  survival  function  for  the  spare  is 


F2(t)  =  e 


-X2t 


t  >  0 


The  survival  function  for  this  system  is 


F(t)  -  F1(t)  +  /  F2(t-s)f1(s)ds 


-A,t  t  -A,(t-s)  -A. s 

=  e  +  f  e  4  A,e  A  ds 


-Ait  \  _“x2t  t 


X1  '  X2 


-(A  -A  )s 

/  (A1-A2)e  1  4  as 


”  A.t  A,  — A«t  —  (A..— A^)t 

-  e  1  +  ■? — =-r —  e  4  (1  -  e  1  4  ) 

A1  "  A2 


For  n  =  3 : 


The  resulting  survival  function  for  the  two  component 
system  will  be  convolved  with  the  survival  function  of  the 
next  component  to  obtain  the  survival  function  for  the  three 
component  system. 

The  result  is 


F(t)  = 


A2A3 


A1A3 


“Alt 

U3-X1)  (A2-X1)  e  +  (X3-X2)  (X1-X2) 


-A2t 


A1A2 


“A3t 


(A3‘A2}  (A3“Al} 


For  all  n: 

Continuing  in  this  manner  one  gets  the  general  solution 
for  the  survival  function  of  this  kind  of  system: 


F  (t ) 


n 


l 

i=l 


n 

XL 

n 


rcf-rrr 


The  two  general  survival  functions  for  systems  built  from 
components  which  all  have  dissimilar  failure  rates  or  all 
have  the  same  failure  rate  are  very  useful  for  computations 
of  survival  functions  of  systems  which  are  formed  from  sub¬ 
systems  of  the  above  two  kinds. 

It  will  be  shown  in  the  following  section  how  to  use 
the  above  results  to  get  a  closed  form  for  the  survival 
function  of  a  system  which  can  not  be  computed  with  only  one 
of  those  general  formulas . 
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3 .  Systems  With  Identical  and  Dissimilar  Components 
Shorthand  notation: 

EXP  [  A^]  +  EXP  [ A x ]  +  EXP[X2J  +  EXP[A2] 

Let  the  random  variables  T  be  distributed  in  the  following 
way: 


Tx  -  EXPUJ,  T2  ~  EXPUJ, 

T3  ~  EXPUj],  T4  -  EXP  [ X 2 ]  . 

There  are  different  approaches  available  for  computing  the 
survival  function  for  this  system. 

One  can  compute  the  survival  function  for  the  first 
two  and  then  the  second  two  components,  using  the  general 
formula  of  Section  III.B.l  and  then  convolve  those  two  sur¬ 
vival  functions;  or  one  can  compute  the  survival  function  of 
the  first  two  components  and  then  add  the  other  components 
one  at  a  time.  The  latter  approach  will  be  used  here  to 
give  a  general  idea  of  the  procedure,  and  to  lead  to  an 
algorithm  for  the  numerical  evaluation  of  the  reliability  of 
systems  consisting  of  so  many  components  that  a  closed  form 
is  no  longer  attainable  or  not  attainable  in  an  efficient 
manner. 

The  first  two  components  (T£  =  t^,  T£  =  T1  +  T2) 


have  as  shown  in  III.B.l  the  survival  function 


( 


_  _  t 

FT,(t)  =  FT,(t)  +  /  F—  (t-s)f,_,(s)ds 

2  A1  0  1 2  rl 

“Xlt 

=  (a11  +  a12t)e  ,  t>0 


with: 


flnd  ci^2  A  ^  i 


Adding  the  next  component  leads  to  +  t2 

and  the  survival  function 


FT^(t) 


FT,(t)  +  /  F-, (t-s) f_, (s)ds 


0  T3 


-A, t  t  -A  (t-s)  _  -A,s 
=  (1  +  A, t) e  +  /  e  2 

A  0 


A^se  1  ds 


which  leads  after  the  integration  to 


ft, (t) 
A3 


(A2  -  2A1A2)  +  (A1A2  -  A^A2)t  -Ajt 
- 5 -  e  + 

(ArA2)2 


(ArA2) 


or  in  the  form  of  equation  (1)  to 


Ft. (t) 
a3 


-A.t  -A2t 

( aii+ai2fc^  e  +  a21e 


Adding  the  last  component  leads  to  T!  =  T,  +  T„  +  +  T,  and 

4  12  3  4 

the  survival  function: 


F«,(t) 

t4 


FT,(t)  +  /  P  , (t-s) f_, (s)ds 


-A.t  -At 

(all+a12t)e  +  a21e 


t  -A2(t-s)  -A.s 

*  J  e  t<al2-allVa12xls)e 


-A2s 

-  ^2a21e  Ids 


Integrating  and  collecting  terms  gives  the  final  survival 
function  for  the  whole  system: 


_  -A.t  -A.t 

FTj(t)  -  (a11+a12t)e  A  +  (a21+a22t)e  ,  t  >_  0 

with 


lll 


a£-3A*A2 
(A2  -  A1) 


‘21 


X1  "  3X1X2 
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A  A2 
A1A2 


(a2-a1) 


22 


X1X2 


(al-a2) 


This  section  demonstrates  that  one  can  get  the  survival 
function  by  adding  one  random  variable  at  a  time  and  thus 
computing  the  new  survival  function.  Each  summation  then 


gives  a  new  or  updated  a_ 


Even  though  the  integration  is 


not  shown  explicitly,  one  can  imagine  that  this  method  of 
convolving  survival  functions  is  very  tedious  and  it  is  very 
easy  to  make  mistakes.  So  the  next  section  introduces  an 
algorithm  to  compute  all  the  a^j ,  which  then  makes  it  easy 
to  evaluate  the  survival  function  at  any  time. 

C.  ALGORITHM  TO  COMPUTE  THE  CONVOLUTION  OF  EXPONENTIAL  LIVES 
The  shorthand  notation  to  be  solved  is 

EXPUJ  +  EXPUJ  +  ...  +  EXPU^  ^  terms) 

+  EXPU21  +  EXPU2]  +  ...  +  EXP[X21  (n2  terms) 

+  . 

+  EXP[Xk]  +  EXP[Xk]  +  ...  +  EXP[Xk)  (nk  terms) 

The  general  solution  for  the  survival  function  is  of  the  form 


F(t) 

-A.t  -A,t 

-  Ax(t)e  A  +  A2(t)e  *  +  ... 

-A.t 

+  Vt)e  K  ,  t  >  0 

with 

Ax(t) 

"  *11  +  *12t  +  *13t2  +  *14t3  + 

(n  -1) 

•**  +  aln1t 

A2(t) 

*  *21  +  *22t  +  *23t2  +  *24t3  + 

(n2_1) 

+  a2n2t 

2  3  ‘V1* 

Ak(t)  -  +  *k2t  +  ak3t  +  ak4fc  +  •••  + 

An  example  shows  how  to  get  the  survival  function  for  the 
following  shorthand  notation: 


EXPl^]  +  EXPU^  +  EXP  [ A x ] 

+  EXPIXj]  +  EXP[X21 

+  exp[x3] 

Survival  Function: 

_  2  -X.t  -X,t  -X^t 

F(t)  »  (a1i+a12t+a13t  )e  +  (a21+a22t)e  *  +  a31e  J  , 

t  >  0. 

An  algorithm  to  compute  all  the  factors  a^,  a12,  . ..,  a31 
and  thus  make  it  possible  to  convolve  all  the  different 
exponential  functions  was  earlier  introduced  by  Gursel  [Ref. 
1] ,  and  is  given  here  in  a  slightly  modified  form.  This 
algorithm  adds  one  exponential  random  variable  in  each  run. 
The  above  example  can  be  solved  in  this  way  by  looping  six 
times  through  this  algorithm. 

Notation  for  the  algorithm: 

K  *  number  of  dissimilar  failure  rates  (3  in  above 
example) ; 

li  -  i-th  type  failure  rate,  i  *  l,...,k; 

nn.  *  number  of  random  variables  having  the  9-th 
1  type  failure  rate; 

1.  »  failure  rate  of  the  ie-th  lambda,  ie  »  1,...,K 

re 

n,  *  current  number  of  random  variables  having  the 
i-th  failure  rate,  n^  =  l,...,nn^? 

a^j  ■  j-th  coefficient  of  the  i-th  polynomial 

ie  -  index  for  current  polynomial  and  currently 
entering  life. 
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Algorithm  to  compute  the  coefficients: 


Initialization:  a.. 


=  0,  for  all  j,k  where  j  =  1,...,K, 
k  =  1, . . . ,nn^ 

=  0,  for  all  i  where  i  *  1,...,K 


Input : 


X^  for  all  i  where  i  =  1,...,] 


nn^  for  all  i  where  i  =  1,.,.,K« 


After  the  first  run:  a11  *  1,  n^  =  1 


n . _  =  n.  +  1,  until  n.  =  nn. 
ie  re  '  re  le 


1.  Update  the  coefficients  a.  ,  for  k  =  2,...,n. 

16  t *  16 

aie,k  *  Xie*aie,m-l/(m-1) '  where  m  -  n.e-j  for 


i  -  0,...,nle-X 


2.  Update  the  coefficients  a.  , 

ie ,  x. 


a  •  i  —  a . 
ie,l  ie 


ni  a. , X  .  ni  (j-l)iX.  a.. 

i  +  l  t  +  l  - i^T1] 

'  i*ie  ui  AieJ  j=2  (X  .  -X  .  ) 3 
i=l  ij^ie  1  ie 


Update  the  other  coefficients  a. .  for  all  i,k  where 

i  <>  ie,  n^  >  0 

for  i  =  1, . . . , K, 
k  =  1,  •  •  •  ,n  . 


ai,nt  -  at,n.Xie/Uie-Xi)  for  a11  1  *  ie' 


(a.  n  X.  ~ma.  m)/(X.  -X.) 

i,n^  le  x,m  ie  i 


for  all  i  where 
i  ?  ie  and  n.  >  1 
for  m  =  n.-j, 
j  =  l,...)(n.-l) 


m 

% 


i 


t  *  | 

.vij 
*  Vi  I 

■A 

!^Vl 


A 
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D.  PROGRAM  FOR  NUMERICAL  EVALUATION  OF  THE  RELIABILITY 
FOR  A  ONE  COMPONENT  SYSTEM  WITH  SPARES 

The  above  algorithm  is  used  to  compute  the  reliability  of 
a  one  component  system  with  similar  and/or  dissimilar  spares 
for  any  given  time.  The  input  to  the  program,  which  is 
written  in  PASCAL  and  listed  in  Appendix  A,  will  be  explained 
here  with  the  same  example  as  above: 

EXP  [  A  x  ]  +  EXPtXj  +  EXPfA^ 

+  expu2j  +  exp[x21 
+  expU3] 

Let's  assume  that  X1  -  0.1,  X2  =  0.2,  X3  ■  0.3  and  we  want 
to  know  the  reliability  for  t  *  0,  10,  20,  . . . ,  90  hours.  The 
prompts  of  the  program  are  typed  in  large  letters. 

"HOW  MANY  DISSIMILAR  LAMBDAS  DO  YOU  WANT  TO  USE?" 


"WRITE  DOWN  ALL  THE  DISSIMILAR  LAMBDAS  WITH  AT  LEAST  ONE 
BLANK  IN  BETWEEN  (INPUT  NUMBERS  <  1  WITH  A  LEADING  ZERO) 

0.1  0.2  0.3 

"WRITE  DOWN  THE  NUMBER  OF  EACH  OF  THE  LAMBDAS  WITH  AT  LEAST 
ONE  BLANK  IN  BETWEEN:” 

3  2  1 

"HOW  MANY  TIME  POINTS  DO  YOU  WANT  TO  BE  COMPUTED?  {<=20) 


"AT  WHAT  TIMES  DO  YOU  WANT  THE  SURVIVAL  FUNCTION  TO  BE 
EVALUATED?  (TIME  SHOULD  HAVE  THE  SAME  DIMENSION  AS  THE 
LAMBDAS  AND  SHOULD  BE  AT  MOST  20  INTEGERS  BETWEEN  0  AND  70. 


THE  UPPER  LIMIT  70  IS  REQUIRED  ONLY,  IF  A  GRAPHICAL 
PRESENTATION  IS  DESIRED  LATER) 

0  10  20  30  40  50  60  70  80  90 

"DO  YOU  WANT  A  SIMULATION  TO  BE  DONE  PARALLEL  TO  THE 
COMPUTATIONS?  (SERIOUS  SIMULATION  NEEDS  A  LOT  OF  TIME!)  HIT 
1  FOR  "YES"  AND  "0"  FOR  NO:" 

1 

"HOW  MANY  NUMBERS  DO  YOU  WANT  TO  USE  AS  BASIS  FOR  THE 
SIMULATION? 

1000 


After  the  computations  the  program  will  produce  the  following 
output : 

FOLLOWING  LAMBDAS  WERE  ENTERED:  0.1  0.2  0.3 
#  OF  EACH  LAMBDA:  321 


Results : 


RELIABILITY 

( 

T 

> 

0 

) 

as 

1.000000 

(SIM 

a 

1.000) 

RELIABILITY 

( 

T 

> 

10 

) 

0.995799 

(SIM 

a 

0.996) 

RELIABILITY 

( 

T 

> 

20 

) 

- 

0.925259 

(SIM 

* 

0.925) 

RELIABILITY 

( 

T 

> 

30 

) 

= 

0.739626 

(SIM 

a 

0.728) 

RELIABILITY 

( 

T 

> 

40 

) 

s 

0.508910 

(SIM 

a 

0.489) 

RELIABILITY 

( 

T 

> 

50 

) 

* 

0.311271 

(SIM 

a 

0.334) 

RELIABILITY 

( 

T 

> 

60 

) 

a 

0.174439 

(SIM 

= 

0.173) 

RELIABILITY 

( 

T 

> 

70 

) 

= 

0.091597 

(SIM 

a 

0.068) 

RELIABILITY 

( 

T 

> 

80 

) 

= 

0.045784 

(SIM 

a 

0.036) 

RELIABILITY 

( 

T 

> 

90 

) 

a 

0.022028 

(SIM 

a 

0.020) 

"DO  YOU  ALSO  WANT  A  GRPAHICAL  PRESENTATION  OF  THE  DATA?  (HIT 
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The  graphical  output  will  not  be  shown  here,  but  will  be 
generated  from  the  program  if  desired. 

Some  notes  on  the  program: 

1.  The  program  is  set  up  to  allow  20  dissimilar  lambdas. 

2.  It  is  necessary  to  type  leading  zeros  in  PASCAL  if  one 
inputs  decimal  numbers. 

3.  The  simulation  generates  n  numbers  (in  this  example 
n  =  1000)  for  every  random  variable,  to  produce  n 
observations  on  the  whole  system. 

4 .  Because  of 

-  the  behavior  of  the  exponential  function, 

-  the  sequence  of  the  lambdas  in  the  input, 

-  the  method  of  adding  always  only  one  random  variable 
at  a  time, 

-  the  possible  closeness  of  the  values  of  the  lambdas 
and 

-  the  possibly  large  number  of  the  dissimilar  lambdas, 
the  program  might  have  to  compute  with  very  large  and 
small  numbers  and  could  produce  erroneous  results. 
Therefore  a  simulation  is  provided  within  the  program 
to  see  if  the  computations  are  close  to  the  simulation. 
Fixed  limits  which  avoid  errors  can't  be  given,  because 
all  the  possible  errors  are  closely  related.  Fewer 
dissimilar  lambdas  for  example  would  allow  that  their 
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values  could  be  closer  together,  or  that  less  care  is 
needed  with  the  sequence  in  the  input,  etc. 


But  after  a  long  test  phase  and  a  comparison  with  results , 
produced  with  the  IBM  3033,  some  hints  can  be  given: 

1.  The  input  sequence  of  the  lambdas  should  go  from  the 
smallest  to  the  largest, 

2.  There  should  be  no  more  than  20  lambdas  altogether, 

3.  The  difference  in  value  of  the  dissimilar  lambdas  should 
be  larger  than  0.05. 
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IV.  RELIABILITY  FOR  A  SYSTEM  OF  ACTIVE  COMPONENTS  AND  SPARES 

A.  MIXTURE  OF  CONVOLUTIONS 

Before  solving  systems  with  active  components  and  spares, 
the  mixture  of  convolutions  will  be  explained  briefly.  Let's 
assume  that  we  want  to  compute  the  reliability  of  a  parallel 
system  with  two  different  components,  like  that  in  Figure  1. 


Figure  1.  Simple  Parallel  System 

The  system  starts  with  a  life  distributed  as  EXP[A^+X2] .  Then 

there  are  two  different  possibilities:  either  component  #1 

fails  with  probability  P-^  =  *l/^l+^2^  and  the  residual  life 

is  distributed  as  EXP[A21«  This  is  the  upper  branch  of 

Figure  2.  Or  component  #2  fails  with  probability 

P2  *  A2/^*l+*2*  and  the  residual  life  is  distributed  as 

EXP  I A  x ]  . 
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EXPlA^+Aj]  + 


A  ^+A2 


exp[a2] 


EXP  f  A  x ] 


Figure  2.  Shorthand  Notation  for  Parallel  System 


Using  the  "MIX "-notation,  the  life  can  be  written  in  the 
shorthand  notation  as: 


EXP[Ax+A2]  +  MIX{P1*EXP[A2] ,P2*EXP[A1] } 


Because  of  the  distributive  law,  which  is  valid  for  these 
mixtures,  one  can  transform  this  equation  to 


MIX{P1*EXP[A21  +  EXP[A1+A21) ,P2*(EXP[A1]  +  EXPIAj+Aj])} 


The  first  and  the  second  expressions  in  parentheses  multiplied 
by  their  respective  probabilities,  is  a  mixture  of  convolu¬ 
tions.  The  following  section  demonstrates  how  those  mixtures 
can  be  solved  numerically  even  if  an  analytical  solution  is 
not  feasible. 


B.  DESCRIPTION  OF  THE  SYSTEM 

This  section  is  concerned  with  computing  the  reliability 
of  a  system  of  components.  The  different  components  might 
be  supported  by  spares,  which  have  the  same  or  a  different 
failure  rate  as  the  original  component.  The  following  example 
will  be  used  to  illustrate  the  idea  of  describing  the  system 
in  the  reliability  shorthand  notation.  The  listing  of  a 
PASCAL  program  for  these  types  of  systems  is  provided  in 
Appendix  B.  The  program  incorporates  the  method  of  program 
#1  to  add  the  components,  one  at  a  time,  to  get  first  the 
convolution  of  the  exponential  distributions  and  then  via 
the  mixture  the  reliability  of  the  whole  system. 

Example: 

The  system  consists  of  two  parallel  components  (#1  and 
#2)  in  series  with  a  third  component  (#3) .  A  spare  is  pro¬ 
vided  for  component  #1.  The  components  have  the  failure  rate 
*1'  *2  an^  ^3  resPect;*-ve^Y*  The  spare  has  the  same  failure 
rate  as  the  component  it  can  replace. 


Figure  3.  Combined  Serial/Parallel  System  with  Spare 
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C.  SHORTHAND  BRANCHING  DIAGRAM 


Xx  /EXP[X1+X2+X31+ 


I-TI^I-/EXP[X2+A3] 


X^+X2+X3 


X^7^EXPtX1+X3] 


exp[x1+x2+x3]+  x 1+x3]+  <- 


Xl+X3 


1  2  a3 


‘ZERO  (6) 


Figure  4 .  Shorthand  Branching  Diagram  for  Serial/ 
Parallel  System  with  Spare 


The  numbers  in  the  numbered  parentheses  in  the  figure 
above  describe  the  different  endpoints  (leaves)  in  this  nota¬ 
tion.  The  connections  between  the  original  system  (root) 
and  the  leaves  will  be  called  branches.  In  the  following: 
the  connection  from  the  root  to  leaf  #1  is  called  branch  #1, 
from  the  root  to  leaf  #2  branch  #2,  etc. 

The  original  system  has  a  life  distributed  as  EXP IX1+X2+X^J 
Three  different  things  could  happen  now.  Either  component 
#1  fails  with  probability  X1/(X1+X  +X3)  or  component  #2  fails 
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with  probability  X2/ (X^X  2+X  3)  or  component  #3  fails  with 
probability  X^/ (X^+X2+X3) .  Let's  continue  in  branch  #1  and 
assume  component  #1  failed  with  the  above  probability .  Now 
the  spare  can  replace  the  failed  component  and  the  system's 
life  is  again  distributed  as  EXP [X^+Xj+X^] .  The  system  is  now 
in  the  same  state  as  at  the  beginning  and  the  same  arguments 
hold.  Let's  assume  new  component  #1  fails  with  probability 
*l^^l+*2+^3^  '  fc^en  remaining  life  is  distributed  as 
EXP[X2+X31 .  If  component  #2  fails  (the  probability  of  this 
failure  is  X2 (X^+X2+X3) ) ,  the  remaining  life  is  distributed 
as  EXP[X^+X3].  The  third  possibility,  component  #3  fails  with 
probability  X 3/ (X^+X2+X3) ,  brings  the  failure  of  the  whole 
system,  denoted  by  the  "ZERO "-notation.  This  says  that  there 
is  a  remaining  life  of  zero.  The  reason  why  this  third  branch 
isn't  simply  omitted  is  that  even  if  there  is  no  remaining 
life,  there  is  a  specified  probability  that  this  event  occurs. 
And  this  has  to  be  taken  into  account  to  compute  the  total 
life  of  the  system. 

The  other  three  branches  of  the  above  figure  are  developed 
analogously. 

D.  THE  SURVIVAL  FUNCTION 

Let's  denote  the  survival  function  of  branch  #1  by  F^t), 
and  the  probability  that  the  events  in  branch  #1  occur  in 
this  special  sequence  as  P^.  Then  the  life  of  branch  #1  is 
distributed  as  P1* (EXP [X1+X2+X3]  +  EXP [X^+X2+X3J  +  EXPIX2+X3J) 


and  the  survival  function  of  this  system  is  the  following 
mixture : 


F ( t)  *  MIX{p1*F1(t) ,p  *F  (t) ,p3*F3(t) ,P4*F4(t) ,p5*F5(t) ,pg*i 


For  the  computation  of  the  coefficients  for  every  single 
survival  function  the  algorithm  of  Section  III.C  can  be  used 
again . 

E.  PROGRAM  FOR  THE  NUMERICAL  EVALUATION  OF  THE  RELIABILITY 
FOR  SYSTEMS  WITH  SEVERAL  ACTIVE  COMPONENTS  AND  SPARES 

This  section  describes  the  program  for  computing  the 

reliability  for  systems  as  described  above.  While  reading 

this  part,  one  should  look  at  Figure  4.  The  program  is 

interactive.  Only  the  input  of  Figure  4  will  be  explained 

here. 

Prompt _ Input _ 

"How  many  branches?"  6  CR  (.Carriage  return) 

1.  Branch  : 


1. 

Element 

• 

• 

2 

CR 

3 

CR 

0 

CR 

2. 

Element 

• 

• 

1 

CR 

1 

CR 

2 

CR 

3 

CR 

0 

CR 

3. 

Element 

• 

• 

1 

CR 

2 

CR 

3 

CR 

0 

CR 

4. 

Element 

• 

• 

1 

CR 

1 

CR 

2 

CR 

3 

CR 

0 

CR 

5. 

Element 

• 

• 

1 

CR 

2 

CR 

3 

CR 

0 

CR 

6. 

Element 

• 

• 

0 

CR 

2. 

Branch  : 

1. 

Element 

• 

• 

1 

CR 

3 

CR 

0 

CR 

2. 

Element 

• 

• 

2 

CR 

1 

CR 

2 

CR 

3 

CR 

0 

CR 

3. 

Element 

• 

• 

1 

CR 

2 

CR 

3 

CR 

0 

CR 

6(t>> 


The  branches  in  the  figure  count  from  top  to  bottom  and  the 
elements  in  each  branch  from  the  leaf  to  the  root.  The  first 


element  in  the  first  branch  is  EXP[A2+A3]  and  one  inputs  only 
the  two  indices  of  the  lambdas.  The  second  element  in  the 
same  branch  is  A^/(X^+X2+^3)  and  one  inputs  first  the  index 
of  the  lambda  in  the  numerator  and  then  the  indices  of  the 
denominator.  After  each  element  input  a  "0".  If  there  is  no 
sixth  element  as  in  branch  #1  to  branch  #5,  no  fourth  element 
as  in  branch  #6,  or  no  element  as  the  ZERO  elements  in 
branches  #3,  #5,  and  #6,  simply  input  a  "0".  Because  of  the 
data  structure  used  in  the  PASCAL  program,  the  system  always 
has  to  be  split  up  into  pieces,  so  that  there  is  only  a  single 
failure  in  each  step.  Because  of  the  algorithm  used,  there 
shouldn't  be  more  than  20  elements  in  each  branch  to  avoid 
errors . 
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V.  SUMMARY 


The  reliability  shorthand  is  a  convenient  way  to  describe 
the  life  of  different  systems.  This  shorthand  then  is  the 
underlying  data  structure  to  compute  numerically  the  relia¬ 
bility  for  systems  for  which  it  is  difficult  or  impossible 
to  get  the  survival  function  in  closed  form. 
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APPENDIX  A 


PROGRAM  CONVOLUTION; 

USES  TRANSCEND , APPLESTUFF ; 

CONST 

ASTERISK-' *' ; 

BLANK  »'  ’ ; 

DASH  ; 

LINE 

TYPE 

ARR  -  PACKED  ARRAY El.. 20 3  OF  REAL; 

ARI  —  PACKED  ARRAY Cl.. 20 3  OF  INTEGER; 

VAR 

As PACKED  ARRAYC1. .20, 1. .203  OF  REAL; 

L . PRO ,  REL  , RELS I M  sARR; 

NINIT,NI, T: ARI ; 

KK, JJ, IE, K, TMAX, SIMMAX, Y, X. COUNT, VAL: INTEGER; 

I , I 1 , J,FACT. NN, NKK, NNI , CHGR, CHSIM: INTEGER; 

IT,D  : INTEGER;  (*  INDEX  FOR  T-ARRAY  *) 

SUM, REL I: REAL; 

PRINTVARs  CHAR; 

<*  L  -  ARRAY  OF  DISSIMILAR  FAILURE  RATES  *) 

(t  NIC  13  «  NUMBER  OF  LIVES  HAVING  THE  ITH  FAILURE  RATE  *) 
(*  AC  I , J3  -  COEFFICIENT  FOR  THE  POLYNOMIAL  *) 

(*  T  -  ARRAY  OF  TIMES  AT  WHICH  THE  RELIABILITY  WILL  BE 
COMPUTED  t) 

PROCEDURE  ATTENTION; 

VAR 

DURATION, Is  INTEGER; 

PITCH: ARRAYC 1. . 103  OF  INTEGER; 

BEGIN 

PITCHC 13s-  15; 

PITCHC23 : —  17; 

PITCHC33  s  —  19; 

PITCHC43  s  —  20; 

PITCHC53:-  22; 

DURATIONS-  30; 

FOR  1 5-1  TO  5  DO 

NOTE (PITCHC I  3, DURATION) ; 

END; 

PROCEDURE  DECIDE2; 

BEGIN 

WRITELN; WRITELN; 

WR I TELN ( ' DO  YOU  WANT  TO  CONTINUE  WITH  THE  SAME  PROGRAM? 
TYPE  "U"  s’); 

WRITELN ('DO  YOU  WANT  TO  COMPUTE  MIXED  CONVOLUTIONS? 

TYPE  "XP2"  s'); 

EXIT  (PI ) ;' 


END; 


PROCEDURE  INPUT; 

BEGIN 

WRITELN; WRITELN; 

G0T0XY(0, 12) ; 

WRITELN (’HOW  MANY  DISSIMILAR  LAMBDAS  DO  YOU  WANT  TO  USE?  ’) 
READLN(K) ; 

WRITELN; 

WRITELN (’WRITE  DOWN  ALL  THE  DISSIMILAR  LAMBDAS  WITH  AT 

LEAST  ONE  BLANK  IN  BETWEEN: ' ) 5 

WRITELN (’ (NUMBERS  <  1  NEED  A  LEADING  ZERO!)’); 

FOR  1**1  TO  K  DO 
READ (LE 13); 

WRITELN; 

WRITE (’WRITE  DOWN  THE  NUMBER  OF  EACH  OF  THE  LAMBDAS  WITH  AT 
LEAST  ONE  BLANK  IN’ ) ; 

WRITELN (’BETWEEN: ’ ) ; 

FOR  1**1  TO  K  DO 
READ (NINITE  ID; 

WRITELN; 

WRITE  (’HOW  MANY  TIME  POINTS  DO  YOU  WANT  TD  BE  COMPUTED"1 
(  <-  20  >  :  ’  )  ; 

READLN(TMAX) ; 

WRITELN (’AT  WHAT  TIME  DO  YOU  WANT  THE  SURVIVAL  FUNCTION  TO 
BE  EVALUATED?’ ) ; 

WRITELN (’ (INTEGER  TIMES  BETWEEN  0-60  IN  THE  SAME  UNITS  AS 
THE  LAMBDAS) ’ ) ; 

FOR  I :*1  TO  TMAX  DO 
READ (TCI  3) ; 

WRITELN; 

WRITELN (’DO  YOU  WANT  A  SIMULATION  TO  BE  DONE  PARALLEL  TO 
THE  COMPUTATIONS?’ ) ; 

WRITELN (’ (SERIOUS  SIMULATION  NEEDS  A  LOT  OF  TIME!)’); 

WRITE (’HIT  "1H  FOR  YES  OR  "0"  FOR  NO  :  ’>; 

RE ADLN ( CHS I M  > ; 

WRITELN; 

IF  CHSIM  *  1  THEN 
BEGIN 

WRITELN (’HOW  MANY  NUMBERS  DO  YOU  WANT  TO  USE  AS  BASIS 
FOR  THE  SIMULATION?’ ) ; 

READLN (SI MMAX ) ; 

END; 

WRITELN; WRITELN;  , 

END; 

PROCEDURE  OUTPUT; 

BEGIN 

WRITE(’ FOLLOWING  ’,K:3,’  LAMBDAS  WERE.  ENTERED  :’); 

FOR  I:*l  TO  K  DO 
WRITE (LCI3:5:2) ; 

WRITELN; 
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WRITE ( *  #  OF  EACH  LAMBDA 

FOR  I:«l  TO  K  DO 
WRITE (NINITCI3:5) ; 

WRITELNj WRITELN; WRITELN; 

•nd; 

PROCEDURE  INITIALIZE; 

BEGIN 

KK:-I; 

WHILE  KK<“K  DO 
BEGIN 

ACKK, 138*1.0; 

NICKK3:*1; 

WHILE  JJ<»K  DO 
BEGIN 

IF  JJOKK  THEN 
BEGIN 

ACKK, 13s -ACKK, 13*LCJJ3/ <LC J J 3-LCKK3 > ; 

END; 

JJ : *JJ+1 ; 

END; 

KKs*KK+l; 

END; 

KKs»K; 

IEs-l; 

END; 

PROCEDURE  SUMMATION; 

BEGIN 

SUM:»0.0; 

FOR  I:«l  TO  K  DO 
BEGIN 

IF  IOIE  THEN 
BEGIN 

SUM:-SUM+ACI, 1 3*LC I  3/ <LC I 3-LC IE3> ; 

IF  NICI3>«2  THEN 
BEGIN 

FACT : *1 ; 

NKK: »NI C 1 I; 

FOR  II: -2  TO  NKK  DO 
BEGIN 

FACT:«FACT*<II-1> ; 

IF  LC I  3  >LC IE3  THEN 

SUM:»SUM+AC I , I I 3*FACT*LC IE3/EXP ( I I *LN ( 
LCI3-LCIE3) > 

ELSE 

BEGIN 

IF  <11  MOD  2>-0  THEN 

SUM:  -SUM+AC  1 ,  1 I 3*FACT*LC IE3/EXP  ( 

I ItLN (ABS (LC I 3—LC IE3) > ) 

ELSE 


\ 
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SUMs-SUM+ACI, II3*FACT*LCIE3/ <-EXP< 

1 1 *LN (ABS (LC I 3—LC IE3  > ) ) > 

END; 

END; 

END; 

END; 

END; 

AC  IE, 1 3 : -SUM+AC IE, 1 3 ; 

END; 

PROCEDURE  SUB1 ; 

BEGIN 

Is-1; 

WHILE  I<-K  DO 
BEGIN 

IF  IOIE  THEN 
BEGIN 

ACI,NICI33:-ACI,NICI3]*LCIE3/ (LEIED-LCI3) ; 

IF  NIC  1 301  THEN 
BEGIN 
J:«l; 

NNIs-NIcn-J; 

WHILE  J<*  <NI C I  3—1 >  DO 
BEGIN 

AC  I « NNI 3 :  =  (LC IE3*AC I , NNI I— NN1 *AC I , NNI  +  1 3 ) / 
(LCIE3-LCI3); 

NNI : «NI C I 3-J; 

END; 

END; 

END; 

Is«I+l; 

END; 

END; 

PROCEDURE  ONEATATIME; 

BEGIN 

IF  NINITCIE3-NICIE3  THEN 
IE:«IE*H 
ELSE 
BEGIN 

NI C IE3: »NI C IE3+1 ; 

Js*0; 

REPEAT 

NN*«NI C IED-J ; 

ACIE,NN3:-LC IED*ACIE,NN-1 1/ (NN-1)  ; 

UNTIL  NN-2; 

SUMMATION; 

SUB1 ; 

END; 

END; 
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PROCEDURE  RAND; 

VAR 

I, II, I I I, COUNT, J: INTEGER; 

R, TSIMsREAL; 

BEGIN 

III:-0; 

COUNT :«0; 

REPEAT 

TSIHj-O.O; 

IIs-1; 

FOR  J:«l  TO  K  DO 
BEGIN 

FOR  I:«l  TO  NINITCJ]  DO 
BEGIN 

R  s “RANDOM/MAX  I NT ; 

TSIH: — TSJM+  <  — (LN (R) /LCII]>  > ; 

END; 

II:«II+1; 

END; 

IF  TSIM>TE IT]  THEN 
COUNT: -COUNT+1; 

IXXa-XXl+lf 
UNTIL  III-SIMMAX; 

RELSIMC  IT] :  -COUNT/SIMMAX  ; 

END; 

PROCEDURE  SUB2; 

BEGIN 

RELI s-O.O; 

IF  TC IT]— O  THEN 
RELI: -1.0 
ELSE 
BEGIN 

FOR  I:-l  TO  K  DO 
BEGIN 

SUM: -0.0; 

NKK: — NI £13; 

FOR  J:-l  TO  NKK  DO 

SUM:«SUM+ACI. J]*EXP<  ( J-l )  *LN (TC IT] )  ) 
PROC I ]:— SUMtEXP (-LC I ]*TC IT] ) ; 

END; 

FOR  I i-l  TO  K  DO 
RELI : — RELI+PROC I ] ; 

END; 

IF  RELI >1.0  THEN 
RELCIT]:— 1.0 
ELSE 

RELCIT]: -RELI; 

IF  CHSIM-1  THEN 
RAND; 


END; 


PROCEDURE  YAXIS; 

BEGIN 

IF  <X“0)  THEN  PR I NT VAR : “L I NE ; 

IF  <<X“0)  AND  <Y*20)>  AND  <TC13“0>  THEN  PR I NTVAR: “ASTERISK; 
IF  <X— 3)  AND  < Y  MOD  2  “  O)  THEN  PRINTVAR: =’ O’ ; 

IF  (X— 3>  AND  <Y«20>  THEN  PRINTVAR: =’ 1 ’ ; 

IF  <X»-2)  AND  <Y  MOD  2*0)  THEN  PRINT VAR:*’ . ’ $ 

IF  (X— I)  AND  <<Y“20>  OR  <Y*0)  )  THEN  PRINTVAR:  *’ 0’ ; 

IF  ( X— — 1 )  AND  <Y*18>  THEN  PRINTVAR: *’ 9’ ; 

IF  (X*— I >  AND  <Y=16>  THEN  PRINTVAR: »’ 8’ ; 

IF  (X— 1)  AND  <Y*14>  THEN  PRINTVAR:  *’ 7’ ; 

IF  (X*-l>  AND  <Y=12>  THEN  PRINTVAR: *’ 6’ ; 

IF  (X— 1)  AND  <Y*10)  THEN  PRINTVAR:  *’ 5’ : 

IF  <X»-1)  AND  ( Y*8)  THEN  PRINTVAR: *’4’ ; 

IF  (X*-l >  AND  <Y-6>  THEN  PRINTVAR: *’3' ; 

IF  <X*-1>  AND  (Y*4)  THEN  PRINTVAR: *’ 2’ ; 

IF  <X*-I)  AND  (Y*2)  THEN  PRINTVAR: *’ 1 ’ ; 

END; 

PROCEDURE  GRAPHICS; 

BEGIN 

COUNT: *1; 

FOR  Y: *20  DOWNTO  1  DO 
BEGIN 

FOR  X: *— 3  TO  76  DO 
BEGIN 

IF  COUNT < =TMAX  THEN 
BEGIN 

VAL: “ROUND <20*RELCC0UNT3> ; 

IF  <VAL*Y)  AND  ( X-TCC0UNT3 )  THEN 
BEGIN 

PR I NT VAR: “ASTER I SK ; 

COUNT:  “COUNTt-1; 

END 

ELSE 

PR I NTVAR : “BLANK ; 

END 

ELSE 

PR I NTVAR : “BLANK ; 

YAX IS" 

WRITE (PRINTVAR) ; 

END; 

END; 

Y:  *0; 

FOR  Xi—3  TO  76  DO 
BEGIN 

IF  (X  MOD  5  “O)  THEN  PR I NTVAR: “LINE; 

IF  COUNT < “TMAX  THEN 
BEGIN 

VAL : -ROUND ( 20*REL  C  COUNT  3 ) ; 

IF  (VAL-Y)  AND  (X-TCC0UNT3)  THEN 
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BEGIN 

PR I NT VAR s -ASTER I SK } 

COUNT: -COUNT +1; 

END 

ELSE 

IF  (X  MOD  5  0  0)  THEN 
PR I NT VAR s -DASH} 

END 

ELSE 

IF  (X  MOD  S  <>  0)  THEN 
PR I NT VAR  s  —DASH ; 

IF  (X— 3>  THEN  PR  I  NT  VAR:  0'  ; 

IF  (X— 2)  THEN  PRINTVAR: ; 

IF  (X— 1)  THEN  PR  I  NT  VAR:  0'  ; 

WRITE (PR I NT VAR) ; 

END; 

WRITE  (’  O  5  10  15  20  25  30  35  40  45 

50  55  60'  )  ; 

WRITELN< ’  TIME  — >’ > ; 

END; 


BEGIN 

INPUT; 

OUTPUT; 

FOR  IT:— 1  TO  TMAX  DO 
BEGIN 

INITIALIZE; 

REPEAT; 

ONEATATIME; 

UNTIL  IE-K+1; 

SUB2; 

IF  CHSIM— 1  THEN 

WRITELNt 'RELIABILITY  <  T  >’ , TE IT3 : 3, ’  )  *  ' , RELE ITU : 8: 6, 
’  (SIM  — ’ , RELSIME IT3 : 5: 3, ’ ) ' ) 

ELSE 

WRITELN(' RELIABILITY  (  T  >' , TE ITD : 3. '  )  =  ' , RELE IT3 : B: 6) 

END; 

WR I TELN ; WR I TELN ; 

ATTENTION; 

WR I TELN ('DO  YOU  ALSO  WANT  A  GRAPHICAL  PRESENTATION  OF 
THE  DATA?' ) ; 

WRITE ('HIT  -1"  FOR  YES  OR  ”0"  FOR  NO  :  ’>; 

READLN (CH6R)  ; 

WRITELN; 

WR1TELN; 

IF  CHGR-1  THEN 
GRAPHICS; 

DECIDE2; 

END. 


APPENDIX  B 


PROGRAM  MI XED_CONVOLUTION; 


USES  TRANSCEND, APPLESTUFF; 

CONST 

ASTERISK-' *’ 5 

BLANK 

DASH 

LINE  »'!’; 

MAX  -10; 

TYPE 

ARR  «  PACKED  ARRAY C 1 . . MAX 3  OF  REAL; 

ARI  -  PACKED  ARRAY E 1 . . MAX  3  OF  INTEGER; 

DARR  -  PACKED  ARRAY  Cl.. MAX , 1 . . MAX  3  OF  REAL; 

VAR 

AR: PACKED  ARRAYCO. . MAX, 0. . MAX , O. . MAX3  OF  INTEGER; 

A: PACKED  ARRAYE1. .MAX, 1. .MAX 3  OF  REAL; 

SUMREL : PACKED  ARRAY  CO.. MAX 1  OF  REAL; 

N: PACKED  ARRAY CO. .MAX 3  OF  INTEGER; 

NINIT, NI ,  T, ZA, INPsARI; 

PROB, RELTOT , LAMB : DARR ; 

L,  PRO,  REL,  Ps  ARR; 

I, J,K,BR, BRANCHES, COUNT, KK, IL,KKK, JJ. IE: INTEGER; 

IT, II,FA,NN,NKK,NNI,CHGR,TMAX,Y,X,VAL: INTEGER; 

SUM, REL 1 , LAM , DEN: REAL ; 

PR I NT V AR , CH : CHAR ; 

(*  L  »  ARRAY  OF  ALL  LAMBDA  t> 

(*  PR0CI3  «  PROBABILITY  OF  EACH  BRANCH  *> 

(*  T  -  ARRAY  OF  ALL  TIMES  AT  WHICH  THE  RELIABILITY  WILL 
BE  COMPUTED  *> 

(*  BRANCHES  *  NUMBER  OF  BRANCHES  IN  SYSTEM  *) 


PROCEDURE  INPUT; 

BEGIN 

WR I TELN ( ' HOW  MANY  DISSIMILAR  LAMBDAS  ARE  IN  THE  SYSTEM?  '); 
READLN(K) ; 

WR I TELN (’WRITE  DOWN  ALL  THE  LAMBDAS  WITH  AT  LEAST  ONE  BLANK 
IN  BETWEEN: * ) ; 

WR I TELN (’ (NUMBERS  <  1  NEED  A  LEADING  ZERO!)'); 

FOR  I :«1  TO  K  DO 
READ(LCII) ; 

WRITELN ( ’ HOW  MANY  TIME  POINTS  DO  YOU  WANT  TO  BE  COMPUTED? 

(  <-  20  >  :  ’); 

READLN ( TMAX ) ; 

WRITELN (’AT  WHAT  TIME  DO  YOU  WANT  THE  SURVIVAL  FUNCTION  TO 
BE  EVALUATED?'); 

WRITELN(' (INTEGER  TIMES  BETWEEN  0-60  IN  THE  SAME  UNITS  AS 
THE  LAMBDAS)’); 

FOR  I i-l  TO  TMAX  DO 
READ(TCID) ; 
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WRITELN; WRITELN; 

END; 

PROCEDURE  ATTENTION; 

VAR 

DURATION, I : INTEGER: 

PITCH* ARRAYC1.. 103  OF  INTEGER; 

BEGIN 

PITCHC13*-  15; 

PITCHC23*-  17; 

PITCHC33**  19; 

PITCHC43*  «  20; 

PITCHC53*-  22; 

DURATION* »  30; 

FOR  Is«l  TO  5  DO 

NOTE (PITCHCI3, DURATION) ; 

END; 

PROCEDURE  OUTPUT 1; 

BEGIN 

WRITE (’ FOLLOWING  * ,K*3, *  LAMBDAS  WERE  ENTERED 
FOR  Is-1  TO  K  DO 
WRITE (LC I  3:5: 2) ; 

WRITELN; 

END; 

PROCEDURE  0UTPUT2; 

BEGIN 

FOR  1**1  TO  BRANCHES  DO 
BEGIN 

FOR  J*-l  TO  MAX  DO 
BEGIN 

FOR  K*»l  TO  MAX  DO 
BEGIN 

IF  (ARC  I,  J,K3-0)  AND  (KOI)  THEN 
ARCI,J,K3:— 1; 

END; 

END; 

END; 

END; 

PROCEDURE  GETDATA; 

BEGIN 

FOR  Is-O  TO  MAX  DO 
FOR  J**0  TO  MAX  DO 
FOR  K*  *0  TO  MAX  DO 
ARCI,J,K3*«- 1; 

NC03S-1; 

WRITELN; 

WRITELN (’HOW  MANY  BRANCHES?’); 

READLN (BRANCHES) ; 

1**0; 


e  u n  «  »aiw  Tv  *i  *r.  • ■r  •*.  a  ••? 


REPEAT 

I:«I+1; 

WRITELN(I:3,  * .  BRANCH:  M; 

Jt-O; 

REPEAT 
J:-J+l  ; 

WRITELN ( Js  3, * .  ELEMENT:  ’>; 

K:-0; 

REPEAT 

Ki-K+I; 

READ(ZACKD); 

ARC  I ,J,K3:»ZACK3; 

NCJ3:«K; 

UNTIL  ARCI,J,K]«0; 

UNTIL  < (ARCI, J-1,NCJ-133«0>  AND  ( ARC  I , J , NC J3 3-0) ) 
AND  (ZAC  1 3*0) ; 

UNTIL  I -BRANCHES; 

0UTPUT2; 

END; 

PROCEDURE  COMPPROB; 

BEGIN 

FOR  I:-l  TO  MAX  DO 
FOR  J:-l  TO  MAX  DO 
PR08CI, J3:«- 1; 

I:-l? 

J:-2; 

WHILE  I <  —BRANCHES  DO 
BEGIN 
DEN: -0; 

K:-2; 

COUNT: -1; 

WHILE  AR C  I ,  J ,  K  3  >0  DO 
BEGIN 

WHILE  ARC  I . J ,  K3  >0  DO 
BEGIN 

DEN: -DEN+LC ARC  I , J , K3 3 ; 

K:-K+l; 

END; 

PROBC I , COUNT] : -LC ARC  I , J , 1 3  3 /DEN; 

COUNT : -COUNT* 1 ; 

J:«J+2; 

K:«2; 

DEN: —0; 

END; 

I:*I+1; 

J:-2; 

END; 

FOR  1:-1  TO  BRANCHES  DO 
BEGIN 
J*-l; 

PCX3:— 1| 
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WHILE  PROBCI, J3>0  DO 
BEGIN 

PCI3*-PCI3*PROBCI,  JD; 
Jt-J+i; 

END; 

END; 

WRITELN; 

END; 

PROCEDURE  COMPLAMB; 

BEGIN 

FOR  Is-1  TO  MAX  DO 
FOR  J*-l  TO  MAX  DO 
LAMBCI,  J3s  —  1; 

u-i; 

J:-l; 

WHILE  I < -BRANCHES  DO 
BEGIN 
LAMs-O; 

Ks»l* 

COUNT: -1; 

WHILE  ARC  I , J , K3  >— O  DO 
BEGIN 

WHILE  ARC  I , J  ,  K3 >0  DO 
BEGIN 

LAM* «LAM*L C ARC  I , J , K 3 3 ; 
K*-K+l; 

END; 

LAMBC I , COUNT  3  s -LAM; 

COUNT s -COUNT* 1; 

J**J+2; 

Ks— 1 ; 

LAM* -O; 

END; 

I*«I*I; 

Js-1; 

END; 

END; 

PROCEDURE  SORT; 

VAR 

TEMP* REAL: 

SORTED: BOOLEAN; 

BEGIN 

REPEAT 

ls-1; 

SORTED* -TRUE; 

WHILE  LCI+13>0  DO 
BEGIN 

IF  LCI3>LCI+13  THEN 
BEGIN 

TEMP*-LCI3; 
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LCI3s-LCI+13; 

LCI+13s-TEMP; 

SORTEOt -FALSE; 

END; 

Is— 1+1; 

END; 

UNTIL  SORTED-TRUE; 

LCI+ns—  1; 

LCI+23s  —  1; 

END; 

PROCEDURE  INITIALIZE; 

BESIN 
KKs  —1 ; 

WHILE  KKC-K  DO 
BEGIN 

ACKK, 13:-1.0; 

NICKK3s-l; 

JJs-1; 

WHILE  JJC-K  DO 
BEGIN 

IF  JJOKK  THEN 
BEGIN 

ACKK, 13S-ACKK, 13tLCJJ3/ (LC J J 3-LCKK3 ) s 
END; 

J J  s  — J  J+l ; 

END; 

KK:-KK+1; 

END; 

KK: — K; 

IEs-I; 

END; 

PROCEDURE  SUMMATION; 

BEGIN 

SUMs  -O.  Of 
FOR  Is -I  TO  K  DO 
BEGIN 

IF  XOIE  THEN 
BEGIN 

SUM: — SU1+AC If13tLCI3/(LCI3-LCIE3>; 

IF  NICI3>-2  THEN 
BEGIN 
FAs-1; 

NKKs— NIEI3; 

FOR  I Is -2  TO  NKK  DO 
BEGIN 

FAs— FAt  < I 1-1 ) ; 

IF  LCI3>LCIE3  THEN 

SUMs -SUM+AC 1 , 1 I 3*FAtLC IE3/EXP  < 1 1 *LN ( 
LCI3-LCIE3) > 

ELSE 
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BEGIN 

IF  (II  MOD  2>-0  THEN 

SUM! -SUM+AC I , I I 3*FA*LC IE3/EXP ( 

1 1 *LN ( ABS (LC I  3— LC IE3 )  ) ) 

ELSE 

SUMs-SUM+ACI, II3*FA*LCIE3/ (-EXP( 

1 1 <LN ( ABS ( L  C  1 3 -L  C  IE  3 ) ) ) ) 

END; 

END; 

END; 

END; 

END; 

AC  IE, 13:-SUM+ACIE, 13; 

END; 

PROCEDURE  SUB1 ; 

BEGIN 

I:«l; 

WHILE  I<«K  DO 
BEGIN 

IF  IOIE  THEN 
BEGIN 

ACI,NICI33:«ACI,NIC233*LCIE3/(LCIE3-LCI3> ; 

IF  NIL  I  301  THEN 
BEGIN 
Js-1; 

NN2s»NICI3-J; 

WHILE  J<* (NIC  I  3-1)  DO 
BEGIN 

AC  I , NNI 3  s  * (LC IE34AC I , NNI 3— NNI *AC I , NNI  +  1 3 ) / 
(LC IE3—LC 13); 

J : «J+1 ; 

NNI s  *NI C I  3— J ; 

END; 

END; 

END; 

I »*I+1 ; 

END; 

END; 

PROCEDURE  ONEATATIME; 

BEGIN 

IF  NIN3TCIE3-NICIE3  THEN 
IEl-IE+1 
ELSE  ' 

BEGIN 

N1 C IE3: «NI C IE3+1 ; 

J;»0; 

REPEAT 

NNl*NI C IE3-J ; 

AC1E,NN3;*LC IE3*ACIE«  NN-1 3/ (NN-1) ; 
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UNTIL  NN-2; 

SUMMATION; 

SUB1; 

END; 

END; 

PROCEDURE  SUB2; 

BE6IN 

RELI s-O.O; 

IF  TIIT3-0  THEN 
RELI: *1.0 
ELSE 
BEGIN 

FOR  X a*l  TO  K  DO 
BEGIN 

SUM: *0.0; 

NKK: *NI C 1 3 ; 

FOR  J: *1  TO  NKK  DO 

SUM:«SUM+ACI, J3tEXP< (J-l > *LN (TC IT3 ) ) 
PROE 1 3 :— SUM*EXP (-LE I 3*TE IT3 )  ; 

END; 

FOR  I:»l  TO  K  DO 
REL I : -REL I +PROC 1 3 ; 

END; 

IF  RELI >1.0  THEN 
RELC IT3:  — 1 . 0 
ELSE 

RELEIT3:-RELI; 

END; 

PROCEDURE  PREPFORLOOP; 

BEGIN 

FOR  I:-l  TO  MAX  DO 
NINITEI3:— 1; 

KKK : * 1 ; 

i*-i; 

IF  LE23<0  THEN 
BEGIN 

NINITE 1 3:— 1 ; 

K:-l ; 

END; 

WHILE  LEI+13>0  DO 
BEGIN 

IF  LC I 3<LE 1+1 3  THEN 
BEGIN 

N1NITEKKK3: *1 ; 

Ks-KKK; 

KKK i -KKK +1 ; 

Ii-Itl; 

IF  LCI+13<0  THEN 
NINITEKKK3:*1 ; 

Ki-KKK; 


uHHru'zasuK.n  »  *  .'•i'*  > '  .». 


A A  A  A  A  .*-  A  .r 


< 


IF 

(X— 1) 

IF 

(X— 1) 

IF 

(X— 1) 

IF 

(  X— - 

-1) 

END; 

AND  <Y-8) 
AND  (Y-6) 
AND  (Y«4> 
AND  (Y«2> 


THEN  PRINTVAR:-' 4*; 
THEN  PRINTVAR:— ’3' ; 
THEN  PRINTVAR:-*2’ ; 
THEN  PRINTVAR:-’!'; 


PROCEDURE  GRAPHICS; 

BEGIN 

COUNT 1*1? 

FOR  Y:-20  DOWNTO  1  DO 
BEGIN 

FOR  X*— 3  TO  76  DO 
BEGIN 

IF  COUNT <-TMAX  THEN 
BEGIN 

VAL: -ROUND (20*RELEC0UNT3 ) ; 

IF  (VAL-Y)  AND  (X-TECOUNT3)  THEN 
BEGIN 

PRINTVAR: -ASTERISK; 

COUNT : —COUNT ♦ 1 ; 

END 

ELSE 

PRINTVAR: -BLANK; 

END 

ELSE 

PR I NTVAR s -BLANK ; 

YAXIS; 

WRITE (PR I NTVAR) ; 

END; 

END; 

Y:  — O; 

FOR  Xs— 3  TO  76  DO 
BEGIN 

IF  (X  MOD  5  -O)  THEN  PRINTVAR: -LINE; 

IF  COUNT <  — TMAX  THEN 
BEGIN 

VAL: -ROUND (20tREL E  COUNT  3  > ; 

IF  (VAL-Y)  AND  (X-TECOUNT3 )  THEN 
BEGIN 

PR I NTVAR : -ASTER I SK ; 

COUNT : -COUNT ♦ 1 ; 

END 

ELSE 

IF  <X  MOD  3  0  0)  THEN 
PR I NTVAR: -DASH; 

END 

ELSE 

IF  (X  MOD  500)  THEN 
PRINTVAR: -DASH; 

IF  (X— *3)  THEN  PRINTVAR: -'  0'  ; 

IF  (X— 2>  THEN  PRINTVAR:-'.*; 

IF  (X— 1)  THEN  PRINTVAR: 0' ; 
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END 

ELSE 

BEGIN 

IIs-I; 

COUNT* *1; 

WHILE  LC  1 1 3-LC  I 1  +  13  DO 
BEGIN 

COUNT* -COUNT* 1* 

NINITCKKK3* -COUNT; 

FOR  IL*«II  TO  9  DO 
LCIL3*-L£IL+13; 

LCHAX3*  — 1| 

K*«KKK; 

END; 

KKK:»KKK+1; 

Is«II*l ; 

IF  (  (  < COUNT >  1 )  AND  (LCI*13<0>  AND  <LCI3>0>>>  THEN 
BEGIN 

NINITCKKK3  * -1 ; 

K:*KKK; 

END; 

END; 

END; 

END; 

PROCEDURE  LOOP; 

BEGIN 

FOR  IT : *1  TO  TMAX  DO 
BEGIN 

INITIALIZE; 

REPEAT 

ONEATATIME; 

UNTIL  IE-K+1; 

SUB2; 

REL  C I T  3 : -REL 1 1 T  3  HP  C BR  3 ; 

RELTOTCBR. IT3*-RELC IT3; 

END; 

END; 

PROCEDURE  YAXIS; 

BEGIN 


IF 

(X-O) 

THEN 

PRINTVAR* -LINE; 

IF 

( <X-0> 

AND 

<Y-20>) 

AND 

(TCI 3-0)  THEN  PR I NT VAR* -ASTERISK; 

IF 

<X— 3) 

AND 

<Y  HOD 

2-0)  THEN  PRINTVAR*  —  *  O' ; 

IF 

(X— 3) 

AND 

(Y-20) 

THEN 

PR I NTVAR  * — * l * ; 

IF 

(X— 2) 

AND 

(Y  HOD 

2-0) 

THEN  PRINTVAR* -' .  '  ; 

IF 

(X— 1) 

AND 

C  (Y-20) 

OR 

(Y-O))  THEN  PRINTVAR* —’ 0’  ; 

IF 

(X— 1) 

AND 

(Y-18) 

THEN 

PRINTVAR*  —' 9” ; 

IF 

<X— 1) 

AND 

(Y-16) 

THEN 

PR I NTVAR  *  —  ’  8 ' ; 

IF 

(X— 1) 

AND 

(Y-14) 

THEN 

PRINTVAR*  —  *  7’ ; 

IF 

(X— 1) 

AND 

(Y-12) 

THEN 

PRINTVAR*-' 6' ; 

IF 

(X— 1) 

AND 

(Y-10) 

THEN 

PRINTVAR* -’S’ ; 
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WRITE (PRINT VAR) ; 

END; 

WRITE  C*  0  5  10  15  20  25  30  35  40  45 

50  55  60*); 

WRITELN (’  TIME  — >’ ) ; 

END; 

PROCEDURE  DECIDE; 

BEQIN 

WRITELN; WRITELN; 

WRITELN (’DO  YOU  WANT  TO  CONTINUE  WITH  MIXED  CONVOLUTIONS? 
TYPE  ,,UM  *’>; 

WRITELN (’OR  DO  YOU  WANT  TD  COMPUTE  CONVOLUTIONS? 

TYPE  "XPl"  :’); 

END; 

BEGIN 

INPUT; 

OUTPUT 1; 

GETDATA; 

OUTPUT? ; 

COMPPROB; 

COMPLAMB; 

FOR  I i-l  TO  MAX  DO 
LCDs- 1; 

FOR  BRs-1  TO  BRANCHES  DO 
BEGIN 
J:*l; 

WHILE  LAMBCBR, J3>-0  DO 
BEGIN 

IF  LAMBCBR, 1 3 >0  THEN 
BEGIN 

LCJ3: -LAMBCBR, J3; 

Js»J+l; 

END 

ELSE 

BEGIN 

LCJ3 s -LAMBCBR, J+l 3; 

J j-J+1; 

END; 

END; 

SORT; 

PREPFORLOOP; 

LOOP; 

END; 

FOR  I*«l  TO  TMAX  DO 
BEGIN 

SUMRELCI3S— 0; 

FOR  BRi-1  TO  BRANCHES  DO 

SUMRELCI3J-RELT0TCBR, I3+SUMRELCI3; 

WRITELN (’RELIABILITY  (  T  >  ’,TCI3s3, ’  >  - 
SUMRELC 1 3 : 8: 6) ; 
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WRITELN; 

ATTENTION; 

WRITELN( 'DO  YOU  ALSO  WANT  A  GRAPHICAL  PRESENTATION 
OF  THE  DATA?’); 

WRITE  ('HIT  "I"  FOR  YES  OR  “O'*  FOR  NO  s  ')? 
READLN(CHGR) ; 

WRITELN; 

WRITELN; 

IF  CHGR-1  THEN 
GRAPHICS; 

DECIDE; 
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